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ABSTRACT 

Context. The presence of dust strongly affects the way we see galaxies and also the chemical abundances we measure in gas. It is 
therefore important to study the chemical evolution of galaxies by taking into account dust evolution. 

Aims. We aim at performing a detailed study of abundance ratios of high redshift objects and their dust properties. We focus on 
Lyman-Break galaxies (LBGs) and Quasar (QSO) hosts and likely progenitors of low- and high-mass present-day elliptical galaxies, 
respectively. 

Methods. We have adopted a chemical evolution model for elliptical galaxies taking account the dust production from low and 
intermediate mass stars, supemovae la, supernovae II, QSOs and both dust destruction and accretion processes. By means of such 
a model we have followed the chemical evolution of ellipticals of different baryonic masses. Our model complies with chemical 
downsizing. 

Results. We made predictions for the abundance ratios versus metallicity trends for models of differing masses that can be used to 
constrain the star formation rate, initial mass function and dust mass in observed galaxies. We predict the existence of a high redshift 
dust mass-stellar mass relationship. We have found a good agreement with the properties of LBGs if we assume that they formed at 
redshift z=2-4. In particular, a non-negligible amount of dust is needed to explain the observed abundance pattern. We studied the 
QSO SDSS Jl 14816, one of the most distant QSO ever observed (z=6.4), and we have been able to reproduce the amount of dust 
measured in this object. The dust is clearly due to the production from supernovae and the most massive AGB stars as well as from 
the grain growth in the interstellar medium. The QSO dust is likely to dominate only in the very central regions of the galaxies and 
during the early development of the galactic wind. 



Key words. ISM: dust, extinction; ISM: abundances; galaxies: ellipticals: chemical abundances, formation and evolution ; individual: 
SDSS J114816.64+525150.3, MS1512-cB58. 
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, 1. Introduction 30 Myr (i.e. the lifetime of a 8 Mq) since the beginning of star 

• • formation, the large observed amount of dust at high redshift 

00 ' Dust modifies the physical and chemical conditions of the inter- ^^^^ ^^^^^^ ^^^^^ p^^j^^p^ ^„ additional dust source. The 

O ; stellar medium (ISM), in particular it prevents stellar radiation g^e II are one obvious stellar site (Maiolino et al. 2004), but it 

O . from penetrating dense clouds, affects the cooling process, and not clear how dust evolution pr oceeds in SN II enviroments 

T^H . depletes elements into grains. Dust absorbs and scatters photons dPwek et aI.Ll2007bHKozasa et aI.Ll2009h . In fact, while submil- 

^ . mosdy at wavelengths < l^m, and reemits them as infrared pho- jj^^^j^^ observations of the Kepler and CasA SN remnants yield 

■ tons. Those processes affect our interpretation of the observa- q. 1-1 M^, of dust (e.g. Morgan et al., 2003, Dunne et al., 2009, 

tional results, for instance those concerning the spectral energy q^^^^ 2009), infrared observations give « 0.1 of 

distribution (SED) and the chemical abundances of the gas. Dust ^^^^ (Helmhamdi et al., 2003, Sugerman, et al. 2006, Kotak et 

is a fundamental component which cannot be neglected espe- 2009; Rho J. et al., 2008, 2009; Meikle et al., 2007; Takaya 

cially when comparing model results to high redshift observa- 2008). A possibility is that such a discrepancy could be 

tions. There are clear evidences that a non-negligible amount of ^^.^^^^ ^^^y. ^ difference in the instrument sensitivity to the 

dust exists at high redshift as indicated by the extinction curves ^-^^^^^^ ^^^^ pj^^^^^ (^^j^ submillimetre, warm for Spitzer). 

for Lyman break galaxies (LBGs) d ShapIey et ali |200lD and by q„ other hand, theoretical calculations may yield up to 1 Mo 

the far infrared (IRVtp millimeter (mm) dust thermal emission ^^^^^ ^j^^ the exact value dependingon the progenitor mass 

of distant quasars (jBeelen et al. , 2006; Bertoldi et al., 2003at chemical composition of the ejectf] (e.g. Dwek, 1998 and 

| WangetaLH2008|) . The main stellar sites for dust production are references therein). However, a large fraction of this dust can 

the cool dense atmospheres of Asymptotic Giant Branch (AGB) destroyed by the reverse shock inside the SN remnant (e.g. 

stars and the supemovae (SNe), and the type of dust formed de- gj^^^j^j ^ Schneider, 2007), with the details depending on the 

pends on the chemical composition of the ejecta. Despite the ^^^^^^ asymmetry in the explosion (Nozawa et al., 2010). 
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fact that AGB stars and SNe la contribute to dust already after 



For instance in Calura et al. (2008) we assume that a 20 Mq star 
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The so-called smoking quasars, namely quasars pr oduced in 
the clouds in the broad line regions (lElvis et al.Ll2002h . might be 
one of the few sources of high redshift dust, but how much dust 
can be formed by the QSO wind is still unknown. Markwick- 
Kemper et al.(2007) determined the composition of dust in the 
broad absorption line (BAL) quasar PG 2112-1-059. They argued 
that the derived crystalline silicate fraction requires high temper- 
atures such as those found in the immediate quasar environment 
in order to counteract the rapid destruct ion from cosmic rays, 
thus lending support to lElvis et al ] (l2002h idea. 

The measurement of metallicity and dust depletion, which 
are products of the (inter-)stellar evolution, provides an impor- 
tant constraint on chemical models. The formation and evolu- 
tion of dust in primordial SNe as well as the du st evolution ef- 
fects o n SED of galaxies have been studied by iNoza wa et al. I 
(120081) and ISchurer et al. (2009J, respectively. Studies of the 
gas, metal and dus t evolution have bee n performed by means of 
analytical models dPwek et al. L 2007ah and simulations (iDwek L 
19981: iMorgariet al.L 120031: llnoueL 120031: 1 Valiante et al"ll2009ah . 
Calura et al. I ( 2008h applied a detailed formulation of dust pro- 
duction and destruction to a detailed chemical evolution model. 
In that paper the model was mostly constrained by means of the 
solar neighbourhood dust determination, and then it was proven 
to be quite flexible, since it can be readily applied to models for 
elliptical and irregular galaxies. The qualitative dust mass evolu- 
tion with time in high mass spheroids was depicted and Calura et 
al. (2008) showed how the dust is an important factor also in the 
very late phases of their lives, namely when the star formation is 
over. In this paper we aim at studying in detail the dust evolution 
in ellip tical galaxies by starting from the model bv lCalura et al. I 
(|2008|) . We update the dust production by changing the dust con- 
densation efficiency in supernova la and II (in order to take into 
account the above mentioned new observational constraints) and 
adding the possibility of QSO dust. We also revise the accretion 
process timescale. We present models of ellipticals of diff'erent 
baryonic masses in the range 10^ - IO'^Mq, and compare them 
with LBGs and with QSO hosts. In particular, after some general 
predictions, we focus on the particular cases: LEG MS 1512- 
cB58(as a proxy for a relatively low mass elliptical) and the host 
galaxy of QSO Jl 148+5251(as a proxy for the most massive 
spheroids). 

The paper is organized as follows. The chemical evolution 
model and the dust model are described in §2 and §3, respec- 
tively. Observational data are briefly summarized in §4 . We pre- 
dict the chemical and dust evolution of elliptical galaxies in dif- 
ferent conditions and apply them to interpret the observational 
data in §5-7. Our conclusions are summarized in §8. 



2. The chemical evolution model for elliptical galaxy 

2. 1 . The chemical evolution code 

The model adopted here is basically the multi-zone model of 
Pipino & Matteucci (2004). The model galaxy is divided into 
spherical shells 0.1 Re/f thick. Most of the results presented be- 
low refer to global galactic properties and are calculated by con- 
sidering all the shells. In practice, quantities like the total mass 
in stars are computed by summing up the single shell contribu- 
tions, whereas the chemical abundances are (gas) mass-weighted 
averages. The galaxy evolves as an open box in which the ini- 
tial gas mass, with primordial chemical composition, rapidly col- 
lapses, on a time scale t, into the potential well of a dark matter 
halo. 



The rapid collapse triggers an intense and rapid star forma- 
tion process, which can be considered as a starburst and lasts 
until a galactic wind, powered by the thermal energy injected 
by stellar winds and SN (la, II) explosions, occurs. At that time, 
the gas thermal energy equates the gas binding energy of and all 
the residual interstellar medium is assumed to be lost. After that 
time, the galaxies evolve passively. 

The equation of chemical evolution for the element / in each 
galactic shell takes the foUowing form: 

dGiit) 



dt 
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Ht - Tm)Qmi(t - Tm)(p{rn)dm 



,dGi(t) 

+( ; )infall , 

dt 



(1) 



where G,(r) - Mgas(t)Xi(t) is the fractional mass of the element 
/ at the time t in the ISM. The quantity X,(f) is defined as the 
abundance by mass of the element ;. By definition = 1. 

The first term on the right side of eq.[T]gives the rate at which 
the element ; is subtracted from ISM by the SF process. The 
variable iff is the star formation rate calculated according to the 
following law: 



(2) 



namely it is assumed to be proportional to the gas mass via a 
constant v which represents the star formati on efficiency. In or - 
der to reproduce the 'inverse wind model' of dMatteucci L 1 1 9941) . 
an earlier version of the now popular "downsizing", we assume v 
as an increasing function of the galactic mass, (see Table[T]i fol- 
lowing the findings of our previous work (Pipino & Matteucci, 
2004). 

The second term is the rate at which each element is restored 
into the ISM by single stars with masses in the range Ml - Mb,,,, 
where M^ is the minimum mass contributing, at a given time f, 
to chemical enrichment and Mb,,, is the minimum binary mass 
allowed for binary systems giving rise to SNIa {3Mq, Greggio 
& Renzini, 1983). The initial mass function (IMF) is 0(ot) oc 
,^^-(1+1.35)^ (Salpeter 1955), and it is normalized to unity in the 
mass interval 0.1 - lOOM©. In particular, Qmi{t - t,„), is a matrix 
which calculates for any star of a given mass m the amount of the 
newly processed and already present element /, which is returned 
to the ISM. The quantity is the lifetime of a star of mass m 
(Padovani & Matteucci, 1993). 

The third term represents the enrichment by SN la for which 
we assume the single degenerate scenario: a C-O white dwarf 
plus a red giant (Whelan & Iben 1973). We refer to Greggio & 
Renzini (1983), Matteucci & Greggio (1986) and Matteucci & 
Recchi (2001) for further details. The predicted SNIa explosion 
rate is constrained to reproduce the present day observed value 
(Mannucci et al., 2008), by fixing the parameter A - 0.09 in eq. 
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The fourth and fifth terms represent the enrichment by single 
massive stars. The initial galactic infall phase enters the equation 
via the sixth term, for which we adopt the formula: 



\ -r — ) infall - -^i,infallLe 

at 



(3) 



where X, intaii describes the chemical composition of the accreted 
gas, assumed to be primordial. C is the normalization constant 
obtained by integrating the infall law over a Hubble time. For 
instance, for a 1O"M0 model, the accretion history is such that 
90% of the initial gas has already been accreted at tg^ (in fact, we 
halt the infall of the gas at the occurrence of the galactic wind). 
Finally, in order to calculate the potential energy of gas, we take 
into account the presence of a dark matter halo. The assumed 
prescriptions are the same as in our previous models (e.g. Pipino 
& Matteucci 2004). 

2.2. Stellar yields 

The yields used in this paper are as follows: 

1. For single low and intermediate mass stars (0.8 < M/Mq < 
8) we make use of the yields by van den Hoek & 
Gronewegen (1997) as a function of metallicity. 

2. For SNIa and SNII we adopt the recently suggested empir- 
ical yields by Frangois et al. (2004). These yields are a re- 
vised version of the Woosley & Weaver (1995, for SNII) 
and Iwamoto et al. (1999, for SNIa) calculations adjusted 
to best fit the chemical abundances in the Milky Way. When 
discussing J 11 48 +5251 we will also present a case in which 
for O and C we have adopted the yields computed assuming 
mass loss in massive stars for Z > as in Maeder (1992). 
In fact, including mass loss in massive stars produces a large 
loss of C and He, thus lowering the O production, and this 
effect is significant only for over-solar metallicity. Recent 
papers (Mc William & al. 2008; Cescutti & al. 2009) have 
shown that these yields are required to reproduce the [O/Mg] 
and [C/O] ratios at high metallicities in the Galactic bulge. 

3. Dust model 

3.1. Stellar dust 



We adopt the dust model of Calura et al. I (l2008h which uses the 
formalism developed by iDwekl (Il998r Let us define XjusiAt) 
as the abundance by mass of the element / at the time t in the 
dust and since G(t) is the ISM fraction at the time f, the quantity 
Gdusr.i - ^dust.i ■ G(t) represents the normalized mass of the el- 
ement / at the time t in the dust. The time evolution of Gdust.i is 
therefore computed as: 



dGdiisijiO 
dt 



HOXdusuit) 



+ r ^|J{t-T,„)S','"Q„^^{t-T„,)(|,{m)dm 
+A I (p{m) 

-0.5 



(1-A) lP(t-T„,)6^'^QUt-T„^)4>(m)dm 

J Mb., 



[ I fin)iffit - T,n2)6iQ,„i{t - T„a)diAdm 

■M„ 
Mb„, 



h(l - A) i/r(f - T,„)6iQ,„i(t - T,„)<f>{m)dm 

Jm,, 

nMu 

Jm^,. 



+ I i//(t - T„)d'/ Q,„i(t - T„,)(p(m)dm 

Gdust.i Gdust.i ^QSO^ I 

+ + Xii//f{t) 

"^destr '^accrJ 



(4) 



where M„, - 8Mo. We refer the reader to Calura et al. (2008) 
for further details and we note that in this equation we neglected 
the wind term which allows us to take into account the fraction 
of dust which is ejected during the galactic wind. The main rea- 
son is that we are focusing on the t < tg^ evolution. The dust 
mass after the wind depends on several assumptions (e.g. fate of 
the dust in the wind stage, wind mass loading, sputtering in the 
hot medium) and we are planning to improve upon the Calura et 
al. (2008) formulation in a future work. Here recall the general 
features of the model as well as the main assumptions. Only the 
main refractory elements, C, O, Mg, Si, S, Ca, Fe, are depleted 
into dust, and we assume that stars can produce two different 
types of grains: i) silicate dust, composed of O, Mg, Si, S, Ca 
and Fe , and i i) carbon dust, composed of C. As suggested by 
iDwek I (Il998h . we consider that the dust producers are low and 
intermediate mass stars, SNIa and SNII. The condensation effi- 
ciencies 6^-^ , 6'" and 6'', for low and intermediate mass stars, 
SNIa and SNII, respectively, are as follows. 

In low and intermediate-mass stars, dust is produced during 
the Asymptotic Giant Branch (AGB) phase (Ferrarotti & Gail 
2006 and references therein). We assume that dust formation 
depends mainly on the composition of the stellar envelopes. 
If Xo and Xc represent the O and C mass fractions in the 
stellar envelopes, respectively, we assume that stars with 
Xo/Xc > 1 are producers of silicate dust, i.e. dust particles 
composed by O, Mg, Si, S, Ca, Fe. On the other hand, C 
rich stars, characterized by Xq/Xc < 1, are producers of 
carbonaceous solids, i.e. carbon dust (Draine 1990). Being 
Mej.iim) oc m X Qmiit) and Md.i(in) the total ejected mass 
and the dust mass formed by a star of initial mass m for the 
element /, respectively, we assume that for stars with Xq/Xc < 1 '■ 

Md.cim) = ■ [M,,;c(m) - 0.75M,y,o(m)] with = 1 

Mdji'n) = 0, for all the other elements. 

For stars with Xq/Xc > 1 in the envelope, we assume: 

Md.cim) = , Md.,{m) = d^'^M.jjim) cc S^'^Q^ 

with = I for Mg, Si, S, Ca, Fe and 

Md,oim) = 16 6^"^ M,jAm)/fi, 

with //, being the mass of the / element in atomic mass units. 
For SNIa we assume: 

Md,i(m) = 6'i"M,jjim) cc 6'^'Q^,,Md,o{m) 
16ZiS'i''M,jj{m)/fXi. 

For SNII, we adopt the same prescriptions as for SNIa. 

The terms Tgccr and Tdestr represent the timescales for the de- 
struction and accretion of dust, respectively. For a given element 
/, the accretion timescale Taccr can be expressed as: 



Taccr.i - Tqj/{1 - fi) 



(5) 
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Table 1. Model parameters and results. 



model 
name 


Ml 


(kpc) 


J' 1 
(Gyr ) 


(Gyr) 


(Gyr) 


L< iyig / 1 € J 


yieiQS 


N 

primary 


ste ai 
dust 


recipe 




QSO 
dust 


(Me) 


(Gyr) 


BEL 
metall. 




M9 


in** 


0.4 


1.2 


0.8 


1 .2 


77^ 

+0.28 


Francois 


no 


off 














M9d 


lU 


0.4 


1.2 


0.8 


1 .2 


+0.28 


Francois 


no 


on 


this work 


0. 1 


\- 

oft 








MIO 


inio 


1 


■J 
J 


U.J 


1 . 1 


_i_n TO 
+u.zy 


Francois 


no 


OIT 














MlOd 


lu 


1 


3 


0.5 


1. 1 


+0.29 


Francois 


no 


on 


this work 


A AO 

0.08 


off 








M3 10 


1 inio 

J ■ lU 


z 


c 

J 


U.J 


T A 
i.U 


_i_A "in 

+U. JU 


Francois 


no 


OIT 














1 i\'K\ 
MJ iUlN 


1 inio 
J ■ lu 


2 


5 


0.5 


1 .0 


+0.30 


Francois 


yes 


off 














Nlj 1 Uu 


J ■ lu 


2 


5 


0.5 


1 .0 


+0.30 


Francois 


no 


on 


this work 


0.05 


off 








Nij 1 UiN a 


J ■ lU 


z 


5 


A ^ 
U.J 


1 A 

1 .u 


_i_n in 

+U. JU 


Francois 


yes 


on 


this work 


U.Uj 


oil 








Nlj 1 UL-d 


J ■ lU 


2 


5 


0.5 


1 .0 


+0.30 


Francois 


no 


on 


Calura 


0.05 


off 








M310MD 


3 ■ 10'° 


2 


5 


0.5 


I.O 


+0.30 


Maeder 




off 














M310MDd 


3- lO'o 


2 


5 


0.5 


I.O 


+0.30 


Maeder 


no 


on 


this work 


0.01 


off 








Mil 


10" 


3 


10 


0.4 


0.9 


+0.33 


Francois 


no 


off 














Mild 


10" 


3 


10 


0.4 


0.9 


+0.33 


Francois 


no 


on 


this work 


0.03 


on 








M12 


10'^ 


10 


20 


0.2 


0.7 


+0.39 


Frangois 


no 


off 














M12d 


10^2 


10 


20 


0.2 


0.7 


+0.39 


Frangois 


no 


on 


this work 


0.01 


on 


2 ■ 10' 


0.049 


mod 


M12MDd 


10'2 


10 


20 


0.2 


0.7 


+0.39 


Maeder 


no 


on 


this work 


0.01 


on 


2 • 10' 


0.049 


mod 


M12dcaseA 


10'^ 


10 


20 


0.2 


0.7 


+0.39 


Francois 


no 


on 


this work 


0.01 


on 


2 • 10' 


0.049 


5x 


M12dcaseB 


10'^ 


10 


20 


0.2 


0.7 


+0.39 


Francois 


no 


on 


this work 


0.01 


on 


4 • 10' 


0.049 


mod 


M12dcaseC 


10'2 


10 


20 


0.2 


0.7 


+0.39 


Frangois 


no 


on 


this work 


0.01 


on 


8 ■ 10' 


0.049 


mod 


M12dcaseD 


10'^ 


10 


20 


0.2 


0.7 


+0.39 


Frangois 


no 


off 




0.01 


on 


1 . 10>* 


0.49 


mod 



The values for f^^,, and [< Mg/Fe >] refer to the central regions of the galaxy. These values are determined only by 
M/,„n (and hence v and r). At each given mass models differ by C, O and N yields - fiducial, Maeder (1992) and primary 
production, respectively - dust (on/off, prescriptions for stellar dust, tq and QSO dust production). The BEL metallicity is 
the metallicity of the gas out of which the QSO dust forms: mod refer to a model self-consistently calculated metallicity, 
5.V is for the case in which the adopted metallicity is 5 times higher than the average mod one. 



where 

Ji = —pr- ■ (6) 

Dust accretion occurs in dense molecular clouds, where volatile 
elements can condensate onto pre-existing grain cores, origi- 
nating a volatile part called mantle (Dwek, 1998, Inoue, 2003). 
Direct evidences for dust accretion come from the observed large 
variations of the depletion levels as a function of density (Savage 
& Sembach 1996) and from the observed infrared emission of 
cold molecular clouds (Flagey et al. 2006), which is character- 
ized by the absence of small grain emission. These features can 
be accounted for by the coagulation of small grains on and into 
larger particles. Indirect evidence for dust accretion comes from 
the estimation of the grain lifetimes, which would be very small 
if no process could allow the grains to recondense and grow 
(McKee 1989, Draine & Salpeter 1979) . 

Dust destruction is primarily due to the propagation of SN 
shock waves in the warm/ionized interstellar medium (McKee 
1989, Jones et al. 1994). Following the suggestions by McKee 
(1989) and D98, for a given element / the destruction timescale 
Tdesir can be expressed as: 

TdestrJ = (eMsNuT^ ' (7) 
'<SN 

Hence the destruction timescale is independent of the dust mass. 
MsNR is the mass of the interstellar gas swept up by the SN rem- 
nant. For this quantity, McKee (1989) suggests a typical value 
of Ms NR ~ 68OOM0 as well as a value for the destruction effi- 
ciency e in a three-phase medium as the present-day local ISM 
are around 0.2, Hence Calura et al. (2008) assumed: 

eMsNR = 0.2 X 68OOM0 = 1360Mo . (8) 

In Calura et al. (2008), Rsn is the total SNe rate at any given 
timestep, including the contributions by both SNIa and SNII. 

Below we discuss the different prescriptions adopted here 
with respect to the original model presented in Calura et al. 



(2008). 

We note here that the modifications have been made through 
changes in the parameters (e.g. S'/, tq,,), whereas the general 
scheme and the above equations hold for both the Calura et al. 
and the present formulations. 
Dust production: 

- We reduced the values for 6'' in order to reproduce recent ob- 
servational constraints (e.g. Kotak et al., 2009, Gomez et al., 
2009). In practice, since the main unknown quantity is the 
efficiency with which newly created dust might be destroyed 
by reverse shocks inside the SN remnant - that we cannot 
model - before it either becomes observable or is ejected in 
to the ISM, we include this uncertainty in the quantity 6'' 
(see also Valiante et al., 2009a, and Zhukovska et al., 2008). 
In particular, we decreased the 6'/ by a factor of 10 with re- 
spect to the Calura et al.'s fiducial case. This implies that a 
typical 2OM0 star now produces nearly O.O8M0 of dust. This 
is 2-3 times higher than the maximum amount of warm dust 
observed in SNII (Rho et al., 2008, 2009), yet lower than the 
most recent estimates of the cold dust mass in the Kepler and 
CasA SN remnants (Gomez et al., 2009, Dunne et al., 2009). 
The dust "yields" as a function of stellar mass for SNell are 
shown in Fig. [1] whereas a comparison between the new 
efficiencies and those in Calura et al. is shown in Table|2] 

- We reduced by a factor of 10 also the values for 6'" and the 
reason is that there is no clear indication that SNe la produce 
dust and certainly they cannot produce more dust than SNe 
II. 

Dust accretion in the ISM: 

We have improved the calculation of the dust growth rate for 
systems with a SFR (and hence a cloud destruction rate) much 
higher than in the MW. 

We assumed that the dust accretion occurs only during the 
starburst epocl£| in dense molecular clouds, where volatile ele- 

^ After the galactic wind occurs and star formation stops there is only 
dust production from long living stars and SNIa. 
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Fig. 1. Dust "yields" from SNe II as functions of initial stellar 
mass adopted in this paper (see description in the text). 



ments can condensate onto pre-existing grain cores. In particu- 
lar, in order to calculate the effective growth timescale we took 
into account the two following facts: i) the growth timescale is 
in the range 10-70 Myr and quickly decreases with the metal- 
Ucity of the ISM (c.f Zhukovska et al. 2008, their Fig. 12), and 
ii) the effective growth timescale depends on the fraction of ISM 
in molecular clouds. As far as the former is concerned, given 
the high SF rates in model ellipticals (see below) the metallic- 
ity quickly increases; therefore, for a sizeable fraction of the 
galaxy evolution, we are in the regime in which tq is shorter 
than the typical survival time of molecular clouds, (20-30 Myr, 
e.g. Krumholz & al. 2006). This implies that the dust growth 
regime is basically determined by the cycling frequency of the 
ISM between clouds rather than by the inner processes occurring 
inside the clouds (Zhukovska et al., 2008), which we will ne- 
glect. Namely, if we consider a fraction of ISM which collapses 
into a cloud which is in turn dispersed by the young stars, there 
is always a net increase in the dust mass. The more frequently 
the single regions of the ISM collapse into clouds, form stars 
and become diffuse ISM again, the faster the gain in dust mass 
by growth. In other words, the only difference between galax- 
ies of different mass/SF history might arise from the fraction of 
ISM in molecular clouds at any time and by the speed at which 
these clouds are (re-)created. In the revised monolithic formation 
scenario, the formation of the most massive galaxies occurs in a 
very intense and short burst, with SFRs exceeding 1000 Mo/yr 
(Fig. 2). Therefore we expect that most of the ISM is in the form 
of clouds for most of the time, during the active star forming 
phase. The high mass galaxy model forms twice as fast (in terms 
of infall timescale, 4 times if the star formation timescale 1/v is 
considered) as the low mass galaxy model. Therefore, relatively 
speaking, we expect that the high mass model will have an effec- 
tive growth timescale which is a factor of 2-4 shorter than in the 
low mass model. In this latter case we adopt the fiducial value 
(to=0.05 Gyr) from Calura et al. (2008), whereas in the former, 
namely the most massive galaxy, we assume to=0.01 Gyr. 
Dust destruction: 

We accounted for the effect of correlated SNe II, which 
strongly suppresses the destruction of the interstellar dust with 
respect to the case in which all SNe could be treated as isolated 
and random explosions. McKee (1989), estimates that in an as- 
sociation with > 10-40 SNell, the dust mass which can be de- 
stroyed amounts to ~ 1.5-2 times the mass of dust destroyed by 
a single isolated explosion, namely the mean efficiency per SNII 
is a factor of 0.05-0.2 the efficiency of an isolated explosion. To 
mimic this effect we made use of an effective SNII rate (equal to 



the 1/10 of the actual SNII rate) when computing the destruction 
rate (Eq. |7]i. The SNIa explosions are not correlated, therefore 
we still used their actual rate for computing the destruction due 
to them. 

Table |2] shows the values adopted for the dust depletion ef- 
ficiencies and both the growth and the destrunction timescale in 
both the old (Dwek et al., 1998, Calura et al., 2008) and in the 
revised (this paper) formulation for the dust evolution. 

3.2. QSO dust 

A non-stellar mechanism for dust production in the early uni- 
verse was proposed by Elvis et al. (2002). These authors showed 
that the physical conditions in the clouds of the Broad Emission 
Line (BEL) regions which undergo cooling and adiabatic ex- 
pansion may become similar to the conditions of AGB stellar 
envelopes, and therefore may provide sites for dust formation. 
The QSO wind can then enhance the pressure, hence boosting 
the dust production. Elvis et al. (2002) estimate that this process 
starts when the BEL clouds are about 3pc from the QSO center. 
Also, they estimate that dust sublimation due to the QSO radia- 
tion should be prevented by the cloud geometry and composition 
even in the most luminous QSOs. 

To take into account this dust production channel w e as- 
sumed a simple model suggested bv lMaiolino et al. I feOOG) with 
the main difference that we use a self-consistent estimate of the 
gas metallicity in the dust formation site. The main ingredients 
are the following: 

1 In order to reproduce the correlation between BH and stellar 
mass, we expect a BH mass M^, ~ O.OOBM/, (Ferrarese & 
Cote, 2007) at the wind time fy,, ■ Since now on, we will re- 
fer to the above "fiducial" set-up, unless stated otherwise. In 
Sec. 6 we will also discuss cases in which the BH-to-stellar 
mass ratio is higher, as it seems to be the case at high red- 
shifts (e.g. Lamastra et al., 2010). 

2 In order to achieve (1) it is enough to let the BH grow at the 
Eddington rate over the time 0-f^,v (Padovani & Matteucci 
1993), starting from a seed mass of ~ 2OOOM0. 

3 We assume that the mass flow rate from QSO wind is i//f{t) = 
0.5 ■ lO-^Mth(t) in units of M^/yr (iProgaet alL|2000i) . 

4 We assume that the QSO dust creation rate scales as the mass 
flow rate in the QSO wind. The factor of proportionality be- 
ing the elemental abundance X, in the gas which is calculated 
at each timestep. Our model predicts average metallicities 
in each shell. Therefore, the abundances in the central region 
of the galaxies are lower than the ones measured in BLRs 
(Nagao et al., 2006b, Juarez et al., 2009). In Sec. 6 we wiU 
show a case in which we estimate the effect of the enhanced 
BLR metallicity in the creation of QSO dust. 

5 We use the same depletion efficiencies as for low- and 
intermediate-mass stars, hence 6'^"° = 6™'. 

Interestingly, several studies (e.g. Maiolino et al. 2004) show 
that the extinction curve in high redshift QSOs is similar to the 
one expected for a medium dominated by SN dust. Indeed, some 
SNell (and possibly "prompt" SNela, see Mannucci et al, 2006) 
can provide the metal seeds out of which the QSO dust can con- 
dense. The effect of the QSO dust amount and composition on 
the spectral properties of high redshift spheroids can possibly 
constrain the above scheme and it will be the topic of a forth- 
coming paper. 

We will study the QSO contribution to the dust only in the 
highest mass galaxy model. The reasons are many: i) the ob- 
servations we compare with are for high redshift very massive 
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Table 2. Dust parameters in the old and the new formulation 



parameter 


Caluia et al. 

1 


new value 
1 


element 


(C/0>1) 








drMg,Si,S,Ca,Ti,Fe 










C 


(C/0<I) 


1 


1 


Mg,Si,S,Ca,Ti,Fe 


6'." 


0.5 


0.05 


C 




0.8 


0.08 


Mg,Si,S,Ca,Ti,Fe^ 


s'.' 


0.5 


0.05 


C 




0.8 


0.08 


Mg,Si,S,Ca,Ti,Fe 


To 


const with gal. mass 


dependent on gal. mass 






true SNII rate 
+true SNIa rate 


effective^ SNll rate 
+true SNla rate 





QSO dust - Elvis et al. (2002) 



" Refer to Sec. 3.1 for the formula for the C dust. 
Refer to Sec. 3.1 for the formula for the O dust. 
A factor of 10 lower than the true rate 



galaxies which host QSOs; ii) the fraction of QSO in z~3 LBGs 
is lower than 3% (Reddy et al., 2008), and iii) there is no indica- 
tion of a QSO in MS 1512-CB58. 



4. Observational Data 

In this section we briefly summarize the observations to which 
we will compare our model predictions in the rest of the paper 



4.1. Lyman break galaxies 

From the discovery of LBGs (see Steidel et al. 1996a,b), only 
a handful of objects whose abundance pattern has been studied 
in great detail. One of the first of such objects for which abun- 
dance measurements were available is MS 1512-cB58, studied 
by Pettini et al. (2002). Owing to its gravitationally lensed na- 
ture, MS 1512-CB58 is one of the brightest known LBGs. This 
object is at z ~ 2.73 and has a luminous mass of ~ lO'^Mf 



a star formation rate iA(f.5/) ~ 40Moyr"' dPettini et al. L l2002h 
and an effective radius of ri ~ 2kpc (Seitz et al. 1998), for a 
Q„, = 0.3, Qa = 0.7, h = 0.70 cosmology. Pettini et al. (2002) 
concluded that the abundances of O, Mg and Si are ~ 2/5 of 
their solar values, whereas [Fe/H]=-1.15, hence underabundant 
by a further factor of > 3 with respect to the other elements. 
This underabundance is probably caused by depletion into dust. 
Pettini et al. (2002) took into account the effect of dust deple- 
tion on Fe-peak elements and suggested that is of the order of a 
factor of two. Matteucci & Pipino (2002) modeled the chemical 
abundances of such a galaxy by taking into account the dust de- 
pletion, as suggested by Pettini et al. (2002) and concluded that 
this galaxy is a small young elliptical undergoing a burst of star 
formation and a galactic wind. We suggested an age of 35 Myr 
for this object. However, that original formulation of the chemi- 
cal evolution model did not take into account the dust evolution 
as we have here. Further confirmation o f the presence of the dust 
came from 3D Lya transfer m odels ( Schaerer & Verhamme L 
I2008h and dust emission models dTakeuchi &IshiiLl2004 

More recently, other lensed LBGs have been observed (e.g., 
Quider et al., 2009, 2010, Dessauges-Zavadsky et al., 2010) but 
the presence of either line emission or intervening systems, as 
well as the different spectral coverage, hampered the measure- 
ment of an extended set of abudance ratios as in MS1512-cB58. 

For the oxygen abundance, we will also use data from the 
AMAZE (Assessing the Mass-Abundance redshift(Z) Evolution, 
Maiolino et al., 2008) and LSD (Lyman-break galaxies Stellar 



populations and Dynamics, Mannucci et al., 2009) programs at 
z~ 3 as well as data from Hainline et al. (2009). 

4.2. Quasars 

The chemical abundances in QSOs are measured either from the 
broad line region (BLR) or from the narrow line region (NLR). 
The BLR traces a region of very small size (of the order of 
parsecs) around the QSO, whereas the NLR extends on sizes 
(of the order of kpc) comparable to the size of the host galaxy. 
The most remarkable facts about abundances in QSOs are: i) 
most of QSOs show oversolar metalU cities, even at very high 
redshift(z > 4) (iMaiohno et al. L l2006h ii) the lack of evolution 
in the QSOs metallicities in a redshift range 2 < z < 4.5 for both 
BLR and NLR (e.g. Nagao et al., 2006a, 2006b, Matsuoka et al., 
2009, Juarez et al., 2009). The most reasonable conclusion from 
this is that the host galaxies of QSO formed at very high redshift 
and evolved very fast due to a very intense star formation which 
produced a lot of metals by means of SNe II. Therefore, over- 
solar metallicities are attained on timescales of hundreds of Myr 
and after that the metallicity cannot increase anymore since the 
host galaxies lose their residual gas by means of galactic winds 
(Matteucci & Padovani, 1993). However, Juarez & al. (2009) 
suggested that the lack of evolution can be due to selection ef- 
fects. 

Re cent data on NLR can be found in iD'Odorico et al. I 
(l2004h . They derived the abundances of C, N and a-elements 
from the NLR of six QSOs. We will also use these data for com- 
parison with our models. 

A non-negligible amount of dust exists in QSO and the dif- 
ferent extinction curve between QSOs and local galaxies are 
challenges for models. 

Here we concentrate o r i the quasar J1148-H5251, wh ich is 
at z ^ 6.4 dBertoldi et al.L l2003bl: llwamuro etaO. |2004 . and 
is one of the most distant QSOs detected (Fan et al. 2003). 
This object is only 375-800 Myr old. Using the FIR luminos- 
ity FfiR ~ IQ^^Lq, the estimated dust mass Mjust ranges from 
2 ■ IQ^Mq to 7 • IO^Mq (e.g. Bertoldi et al., 2003a, Robson et al. 
2004; Carilli et al. 2004; Beelen et al. 2006). 

Several studies (Barth et al. 2003; Maiolino et al. 2005; 
Becker et al. 2006) indicate a near solar metallicity in the quasar 
host. SDSS J1148 shows a bolometric luminosity of Lb„i = 
IO'^Lq powered by accretion onto a supermassive black hole 
(SMBH) of mass 1-5 ■IO'^Mq (Willott et al., 2003, Barth et al. 
2003). The estimated SFR of the host galaxy is ~ 3 ■ 10^ M^yr^ 
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(Bertoldi et al. 2003b; Carilli et al. 2004). Several other high 
redshift quasars are now detected at submillimeter wavelengths 
(Priddey et al. 2003, Wang et al., 2008 and references therein), 
which indicate the presence of dust grains in excess of lO^M©. 

5. Overall results 

The main parameters of the model are shown in Table [1] In 
particular, the model mass is uniquely associated to effective ra- 
dius Reff, Star formation efficiency v, infall timescale r. These 
quantities set the predicted time for the occurrence of a wind 
and the average [Mg/Fe] in the stars. These last two values are 
for the central (0 - Q.IR,,//), because observations typically con- 
strain [Mg/Fe] in such a region. The adopted empirical relation 
between v, t and the galaxy mass is motivated by our previous 
works (e.g. Pipino & Matteucci, 2004). For each model mass, we 
will then consider cases where the C and O yields are taken from 
Maeder (1992), the N yields include primary production in mas- 
sive stars, as well as different dust recipes. These prescriptions, 
however, do not affect the galactic wind timescale and the SFR. 
Here we briefly mention some general properties of the models, 
such as the star formation histories (which comply with down- 
sizing, e.g. Thomas et al., 2005), the dust masses and the abun- 
dance ratio evolution. In Pipino & Matteucci (2004) the mean 
formation redshift of the model galaxies was constrained to the 
observed ages (inferred from optical line-strength indices) and 
the mean colours at a given mass. For the purpose of this paper, 
we leave the formation redshift free in order to find the best fit 
to the single objects that we try to reproduce. 

5.1. Star formation rates 

Fig. |2] shows the predicted SFR for our model galaxies summed 
over the shells and truncated at the average fg„.. It is worth not- 
ing that the average value of f^jt is lower than that calculated for 
the inner regions and reported in table 1 . The reason is that, in 
the adopted framework, elliptical galaxies form outside-in (e.g. 
Martinelli et al., 1998, Pipino et al., 2008), hence a wind devel- 
opes first in the outskirts, whereas the innermost regions are still 
experiencing star formation (and QSO growth). The star forma- 
tion efficiencies have been chosen to reproduce the observations 
- including those in local ellipticals - and are not considered as 
free parameters in the remainder of the paper As discussed also 
in Matteucci & Pipino (2002) and it is also clear from Fig.|2] as 
far as the SF histories are concerned, a model with mass in the 
range 10'" - 10"Mq is the only one which must be compared 
to MS1512-cB58, whereas the most massive galaxy should be 
compared to the host of the QSO J1148+5251. In particular, 
the star formation rate evolution in Fig.|2]shows the downsizing, 
namely that massive galaxies have more efficient star formation 
with a shorter duration than less massive galaxies. 

5.2. Dust masses 

The evolution of the dust mass and of the dust-to-gas ratio are 
shown in Figs. [3]and|4l respectively. A dust mass-galaxy mass 
relation, namely that dust content is higher and increases faster 
in more massive galaxies, is evident. A linear regression analy- 
sis yields the following relations 

Log{Maus,,peaklM^) = O.mLogMUM^) - 2.75 (9) 
when the dust mass at the peak is considered, whereas 
LogiMdusumlM^) = 1 .07(LogMn,JMo) - 4.01 (10) 
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Fig. 2. The predicted star formation rate as function of galac- 
tic age for the model galaxies. Note that the setup for the dust 
and adopted yields do not modify the predicted star forma- 
tion rate. Therefore each line in the plot is representative of all 
models with a given mass. The horizontal lines correspond to 
a SF R=40MHyr~' which is the estimated value for MS 1512- 
cB58 (lPettini et al. Ll2002h . and to SFR=3 x lO^Mojr which is 
the value inferred for Jl 148-1-5251, respectively. 
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Fig. 3. The evolution of the dust mass in the galactic models. 



after 0.5Gyr of evolution. Such a relation stems from the SFR- 
mass relation in Fig.|2]and it is likely to present a redshift evolu- 
tion similar to the observed one in the mass-metallicity relation 
(Maiolino et al., 2008). This is a clear and novel prediction of our 
set of models. Future observations can provide a confirmation 
of this picture. For the 3 ■ 1O'"M0 mass case we also show the 
dust mass which would be obtained by using the prescriptions 
by Calura et al. (2008, thin solid line) compared to a model with 
same mass, but the new prescriptions for the dust (dotted line). 
This latter model will be considered the fiducial case to study 
the galaxy MS 1512-cB58. The Calura et al. prescriptions yield 
a faster and earlier increase in the dust mass, due to the larger 
d'' . At later times, however, the effect of the stronger destruction 
term hampers further increase in the dust mass as opposed to 
the model with the new prescriptions. The implications of such 
a different behavior will be further discussed below by means of 
the chemical abundance ratios. 
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Fig. 4. The evolution of the dust to gas ratio in the galactic mod- 
els. 



5.3. Abundance ratios in tine gas 

The different timescales of the enrichment processes for a - 
elements (massive star ending as SNe II), and iron (dominated 
by intermediate mass stars in binary systems ending as SNe la) 
imply an enrichment delay of iron and this explains the evolu- 
tion of all curves in Fig.|5](upper panel). If we then compare the 
various models in Fig. |5] (upper panel) we can see the effect of 
downsizing very well. Due to the relation between star formation 
efficiency v, infall time scale r and the total luminous mass M/„,„ 
mentioned in §2 , we predict that the [Mg/Fe] ratio in the gas out 
of which stars form is higher for more massive galaxies. This 
ensures the fit of the [ff/Fe] versus mass for ellipticals, as shown 
by Pipino & Matteucci (2004). Moreover (Fig. |5] lower panel), 
the increase in the [X/H] with time is much faster and saturates 
at higher (over-solar) values in the more massive galaxies than in 
low mass ones. This ensures that, despite more massive galaxies 
have less time to convert gas into stars, they do it more efficiently 
than lower mass systems, thus reproducing the mass-metallicity 
relation that we observe in the stars of local ellipticals, namely 
a increase of a factor of ~3 in the metallicity in the mass range 
10^ - IO'^Mq. Until recently, only the mass-stellar metallicity 
relation (e.g. Thomas et al., 2005) could be measured, namely a 
relation that is inferred from the integrated light of local ellipti- 
cals. While important in constraining the galaxy formation pro- 
cess, such a measure reflects an average metallicity and, hence, 
does not allow to capture the details of the star formation histo- 
ries as opposed to the study of the, e.g., O/H ratio versus time. 
Very recent observations (e.g. Maiolino et al., 2008), provided 
the opportunity to directly compare model predictions with gas 
abundances of z~3 star forming galaxies. Indeed Calura et al., 
(2009) find that the observed mass-gas metallicity relation at 
such high redshift can be explained only by models of ellipti- 
cal galaxy progenitors as the ones studied here. Here we do not 
repeat the analysis. Moreover, we show all the models as if they 
were forming at the same time. In principle, both panels of Fig.|5] 
can be used as a diagnostic plot for high redshift galaxies: the 
measurement of the chemical abundance of one element which 
is not easily dust depleted (like O) and the abundance ratio of 
two elements which have very similar condensation efficiencies 
(like Mg and Fe) can immediately tell the age and the mass of 
the galaxy as well as the past SF history, hence the likely mor- 
phology at z=0. 
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Fig. 5. The predicted [Mg/Fe] -[Fe/H] relation for models with 
different masses (upper panel) and the evolution of the metallic- 
ity in terms of the O abundance (lower panel). 



Finally, we show the [N/Fe] (Fig.|6]) and [N/O] (Fig.|7]i ratios 
as a function of the metallicity for galaxies of differing masses. 
The predictions are made for models with (lower panels) and 
without dust (upper panels). Here we anticipate from the next 
section the observed values in the galaxy MS 1512-cB58 , since 
they represent a valu able reference point t o guide the eye. These 
data are taken f rom Teplitz et al. j2000h (Red Cross, T data, 
emission lines), iPettini et al. ( 2002l) (Black Square,P data, ab- 
sorption lines). Only the Pettini et al. data have been used as 
constraint in our previous works. Here we assume that the two 
datasets limit the observed region and the [N/Fe] from Pettini 
et al. as a lower limit which must be necessarily reproduced by 
the models. From a joint analysis of Figs.|6]and|7]we infer that 
only models with mass lower than 3 • lO^^M© predict the [N/Fe] 
within the observed range at [Fe/H]=-1.15. On the other hand, 
only models with mass larger than 3 ■ lO'^M© predict values for 
the [N/O] ratio consistent with observations. Therefore, in agree- 
ment with what has been inferred from the SFR analysis (Fig. 2), 
we will focus on a 3 • 10 "'M© model to study the galaxy MS 
1512-CB58. For this model, the predictions based on the Calura 
et al. (2008) prescriptions are also shown. Since N is not de- 
pleted into dust and O has a very small depletion, both the new 
and the Calura et al. prescriptions lead to basically identical re- 
sults in Fig. I2I On the other hand, the change is more evident 
in Fig. |6] As seen above, in fact, the Calura et al. prescriptions 
give more dust at earlier times with the net effect of increasing 
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the predicted [N/Fe] ratio by 0.3 dex. At [Fe/H]~ -1.15, nearly 
50% of the Fe is locked into dust following the Calura et al. pre- 
scriptions, whereas the fraction is only 12% with the new ones. 

Here it is worth mentioning that [Fe/H] = -1.15 is attained 
at 0.05 Gyr in model Mil, and at 0.15 Gyr in the MIO case. 
In other words, at a fixed metallicity [Fe/H] - -1.15, the model 
Ml 1 has [N/H]< -2, whereas the model MIO has [N/H] ~ -0.8, 
hence a much higher [N/Fe]. This happens because the smaller 
galaxy takes a longer time to attain the same metallicity. A sim- 
ilar explanation holds for the [N/0]-[0/H] diagram. It can be 
shown that this is a consequence of assuming v increasing with 
mass (i.e. the downsizing), whereas, if one assumes that v de- 
creases with mass as in the original monolithic collapse model 
by Matteucci & Pipino (2002), more massive galaxies will ex- 
hibit systematically higher [N/O] values than low mass ones at 
a fixed [O/H] (see, e.g.. Figs. 2 and 4 in Matteucci & Pipino, 
2002). Such a difference in the mass dependence of the pre- 
dicted abundance ratios trends with metallicity is a non trivial 
effect of downsizing that, to our knowledge, has not been shown 
elsewhere and still awaits confirmation by observations. 

We remind the reader should that the uncertainties in the stel- 
lar yields (e.g. Romano et al., 2010) affect the galaxy properties, 
e.g. the age, that we infer when comparing models to observa- 
tions and indeed it is difficult to assign a quantitative measure 
(i.e. an error) to this effect. For instance, since low mass models 
have a subsolar gas metallicity for most of their evolution, we 
expect a negligible impact of the mass-loss dependent O yield 
which seem to be required to explain the abundance pattern in 
metal rich systems as the Galactic bulge (McWilliam & al. 2008; 
Cescutti & al. 2009). Indeed, as it can be seen from Figs.|5]and|6] 
a 3 • 1O"^M0 mass model which differs from the fiducial case only 
for the C and O yields, modified according to Maeder (1992), 
does not predict any significant departure in the abundance ra- 
tios in the metallicity regime which is relevant for our study. We 
do see, however, a ~0.2 dex lower O abundance (lower panel of 
Fig. |5]l with respect to the fiducial case at times later than ~0.6 
Gyr. This happens when the metallicity exceeds the solar value 
(~8.7 - Asplund et al., 2005 - in the units used in the figure) 
and the strong mass loss in the Maeder (1992) yields does not 
allow the O abundance to increase futher We will see that this 
is not the case for the most massive and metal rich objects in 
Sec. 7. In the next section we show some examples of modifica- 
tions to N yields. Here we note that not only the yields, but also 
the prescriptions for binary stars which explode as SNIa may af- 
fect our results. We checked that, given the relatively short time 
needed to attain [Fe/H] — 1.2 in the models, most of the Fe has 
been produced by SNII and changes in the SNIa prescriptions 
do not really affect the predicted abundance ratios. This can be 
easily verified by looking at the upper panel of Fig. |5] all the 
models are still in the plateau phase. Therefore the results in 
this age/metallicity regime that we will present in Sec. 6 are ro- 
bust against changes in the SNIa prescriptions. However, we ex- 
pect more relevant effects at later ages/larger [Fe/H]. Adopting 
recent observationally motivated prescriptions for both prompt 
and late SNIa as in Mannucci et al. (2006), the changes in the 
abundance ratio evolution as a function of [Fe/H] are minimal 
(c.f. Fig. 8 in Matteucci et al., 2006). If instead, one assumes 
the Double-Degenerate scenario (Iben & Tutukov, 1984) as op- 
posed to the Single-Degenerate one adopted here, recent work 
by Greggio (2005) and Valiante et al. (2009b) show that the ex- 
pected number of SNIa event in the former case, and hence the 
Fe produced, is much lower than in the latter Note that, in or- 
der to make a meaningful comparison, all the different channels 
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Fig. 6. The predicted [N/Fe] -[Fe/H] relation for models with dif- 
ferent masses (with and without dust; lower and u pper panel, re- 
spectiv ely). Data for MS 15 12-cB58 are taken fr om lTephtz et a"i~l 
(l2000h (Red Cross, T data). |Pettini et al. I (l2002h (Black Square,P 
data). 

must be constrained to reproduce the observed present-day SNIa 
rate. 

6. The low mass case: comparison with LBGs 

In this section we focus on our low mass model in comparison 
to the class of the LBG galaxies. As noted in the Introduction, 
MS1512-cB58 is still the best candidate for having an exte- 
sive sample of well measured abundance ratios. Therefore, this 
galaxy will is the best case to test our model. 

Absolute abundance evolution usually depends on all the 
chemical evolution model assumptions, whereas the abundance 
ratio evolution only depends on stellar yields, stellar lifetimes, 
and IMF. Therefore, abundance ratios can be used as cosmic 
clocks if they involve two elements enriching the ISM on quite 
different timescales, such as [a/Fe] and [N/a] ratios. These 
ratios, when plotted as functions of metallicity tracers such as 
[Fe/H], allow us to clarify the particular star formation history 
(see Matteucci 2001). 

6. 1 . [N/O] and [N/Fe] ratio 

As shown by Matteucci & Pipino (2002), it is useful to study 
the abundance ratio of a refractory element (like Fe) to an unde- 
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[0/H] 



Fig. 7. The predicted [N/0]-[0/H] relation for models with 
different masses (with and without dust: lower and up- 
per p anel, respecti v ely). Data for MS 1512-cB 58 are taken 
from iTeplitz et al. I (l2000l) (Red Cross, T data), iPettini et al.1 
(l2002l) (Black Square.P data). 




-1.2 -1 

[Fe/H] 

Fig, 8, The predictions of [N/Fe] vs [Fe/H] ratios for a 3 ■ lO'" 
model galaxy. Cases with either primary+secondary or sec- 
ondary only N production in massive stars are shown along 
with models with a nd wi thout dust depletion. D ata are taken 
from iTepUtz et al. I (|2000|) (Red Cross, T data), iPettini et al.1 
(l2002l) (Black Square.P data). 



Fig. 9. The [N/O] vs [O/H] trend. Cases with either pri- 
mary-Hsecondary or secondary only N production are shown 
along with models with a nd without dust depleti on. Data are 
from ITeplitz et al. I (l2000l) (Red Cross, T data), IPettini et al.1 
(l2002l) (Black Square. P data). 



pleted one (like N) in order to measure the amount of dust deple- 
tion for the Fe-peak elements. In Fig. [8] we show the predictions 
from the 3 ■ 10"^ galaxy model with different prescriptions re- 
garding the dust and the N primary production. Fig.[8]shows that 
in order to predict a value for [N/Fe] in the observed range, dust 
depletion is needed, although such a conclusion suffers from a 
degeneracy with the N production. Nitrogen is generally consid- 
ered mainly a secondary element, in the sense that N needs the 
C and O, originally present in the star, to be created during the 
CNO burning cycle. However, N can also be produced as a pri- 
mary element in some special situations, such as in AGB stars 
(Renzini & Voli, 1981; Van den Hoeck & Groenewegen, 1997; 
Marigo 2001) and massive rotating stars (Meynet & Maeder, 
2002). In fact, the C and O going to form N are not the orig- 
inal ones present in the gas out of which the star formed, but 
they have been synthesized by the star itself. In the two differ- 
ent cases N production changes significantly, especially at low 
metallicities. From a look at Figure 9 we can derive the follow- 
ing considerations: it is difficult to say whether the N is primary 
or secondary in massive stars. First of all because Fe is depleted 
into dust and second because of the discrepancy between the N 
measured from emission and absorption lines. On the other hand. 
Figure 10 can give more indications. In particular, it seems that 
assuming primary N from massive stars overproduces N with re- 
spect to O in the gas by more than a factor of three. On the 
other hand, the model which best reproduces the [N/Fe] value 
from Pettini et al. (2002), namely the fiducial model with dust 
and without primary N, is offset by Pettini et al.'s [N/O] ratio by 
0.8 dex. Here we add that at [Fe/H]=-1.15 the model has an age 
of 100 Myr, therefore we confirm that MS 1512-cB58 must be 
young. 

6.2. a - elements to iron ratios 

Fig. [To] and Fig. [TTI show the predicted [Mg/Fe] and [Si/Fe] vs. 
[Fe/H] compared to observational data for the LBG MS 1512- 
cB58 (Pettini et al. 2002). Our predicted [Mg/Fe] and [Si/Fe] ra- 
tios are still lower than observed in MS 15 12-cB58. On the other 
hand, the predicted values are higher by ~ 0.2 dex than predicted 
by Matteucci & Pipino (2002), because the yields adopted in 
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Fig. 10. The [Mg/Fe] vs [Fe/H] ratios for the 3 ■ IO^'Mq galaxy Fig. 11. The[Si/Fe] vs [Fe/H] ratios for the 3 ■ 10'"Mo galaxy 



model. Data from Pettini et al. (2002). 



model. Data from Pettini et al. (2002). 



this work consider a sligthly higher Mg production from mas- 
sive stars with respect to the Woosley & Weaver (1995) yields 
(namely the yields adopted by Matteucci & Pipino, 2002). Note 
that in the so-called "horseshoe" LBG (Quider et al., 2009), the 
Si/Fe ratio is very close to the one observed in cB58, thus con- 
firming the discrepancy between theory and observations. We 
do not believe that the discrepancy could indicate a problem re- 
lated to too low Mg and Si yields in the adopted nucleosynthe- 
sis, since a [Mg/Fe] > 0.6 dex in the gas for a large fraction of 
the galaxy evolution would imply [<Mg/Fe>sM,s] >0.3-0.4 dex, 
which exceeds the observed value observed in present day low 
mass ellipticals. 

When the metallicity is measured from the photospheric ab- 
sorption features of the stellar light (Rix et al., 2004), instead, the 
[Fe/H] is 1.5-4 times higher than those reported by Pettini et al. 
(2002), hence there is the possibility that this disagreement orig- 
inates from the dust, i.e. a possibility could be that in MS 1512 
-cB58 Si and Mg are not as depleted as much as Fe, at variance 
with our prescriptions. For instance, if we leave S'^j^ = S'^'.-O.Odi 
as prescribed by the new model presented in this paper and in- 
crease 6"^ to 0.5, we predict [Mg/Fe]=0.75 and [Si/Fe]=0.6 at 
[Fe/H]=-1.15, hence formally in agreement with the observed 
value, within the error. This, however, cannot be justified on the- 
oretical grounds since we assume the same efficiencies for Mg, 
Si and Fe because they have nearly the same condensation tem- 
perature (see Calura et al.2008 for details). However, such a 
differential depletion effect, which is known to affect the abun- 
dance pattern of the cold gas in Damped Lyman Alpha systems 
(Vladilo 2002, Calura et al. 2003), could partially alleviate the 
discrepancies between models and observations. The difference 
between predicted and observed values is not so evident in the 
Si/C abundance in QSOs (see below). 

6.3. The 0/H ratio and the age of LBGs 

Finally, in Fig. [12] we show the predicted age-oxygen relation 
for Lyman-break galaxies, compared with the recent data from 
MaioUno et al. (2009, AMAZE), Mannucci et al. (2009, LSD), 
Hainline et al. (2009) and Teplitz et al. (2000). Here we obtain 
a good fit to the data if we assume a galaxy formation redshift 
for these galaxies in the redshift interval ~2.3-4. Note that all 
the galaxies have masses around a few 10"^ M©, therefore the 
inferred spread in the formation epoch can be understood as the 
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Fig. 12. [O/H] ratio as a function of redshift for a 3 • 
10'"Mq galaxy and different case formation redshifts. Data are 
from Maiolino et a l (200 9, AMAZE), Mannucci et al. (2009, 
LSDi lPettini et al. I (l2002h (MS 1512-cB58). Clone and Cosmic 
Horseshoe are from Hainline et al. (2009). 

redshift range in which low mass ellipticals form. In particular, 
the curve which gives the best agreement with MS 1512-cB58 
data in Fig. fT2l requires a formation redshift z~3, that is, in the 
standard LCDM cosmology, an age of at least 0.2 Gyr (we recall 
that MS 1512-cB58is at z=2.73). A similar conclusion can be 
reached from the analysis of Fig. 2, where the best agreement 
with the observed SFR is obtained at ages larger than 0.2 Gyr 
However, given the uncertainties in the observational values and 
the flat trend with time in the predicted one, such a constraint 
is rather weak. Finally, we showed above that the time at which 
our model reaches the metallicity [Fe/H] =-1.1 5 is 0.1 Gyr, again 
hinting for a very young age for the galaxy MS 1512-cB58 . 

6.4. Abundance ratios in the dust 

Finally, in Fig.[T3]we show the predicted abundance ratios in the 
dust for the lowest (M9d) and highest (Ml 2d) mass cases com- 
pared to an intermediate mass one (M310d). The high [Mg/Fe] 
and [Si/Fe] are easily explained by the fact that Mg, Si and Fe 
have similar condensation efficiencies, therefore the abundance 
ratios in the dust mirror those in the gas. On the other hand O 
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Fig. 13. Predicted abundance ratios in the dust as function of 
[Fe/H] for a galaxy of 3 ■ lO'^Mg (Model M310d). 



- the dominant metal in the gas - is underabundant with respect 
to Fe in the dust. In passing, we note that the trend in the dust 
[Mg/Fe], [Si/Fe] and [O/Fe] ratios versus [Fe/H] does not sig- 
nificantly with the galaxy mass. We stress that these are actual 
predictions of the models which await observational confirma- 
tion and, if so, will help in understanding the abundance pattern 
in high redshift galaxies. 



7. The high mass case: comparison to a QSO host 

In this section we apply our model to the challenging case of 
the Jl 148+5251 host galaxy. In the first place, we note that our 
model predicts a stellar mass of ~ 1O"M0 and a gas mass of 
~ XO^^Mq within 2.5 kpc (i.e. 1/4 Reff) in agreement with the 
observations within the same aperture (e.g. Wang et al., 2010). 
The predicted BH mass at the onset of the galactic wind is ~ 
2 ■ 10"^ Mq, as required by observations (e.g. Barth et al., 2003). 
In the upper panel of Fig.[T4]we show the predicted dust mass as 
a function of galactic age and compare it with the observed dust 
mass (shaded area). The model, which has been calibrated on the 
LBGs and present-day properties of ellipticals, does reproduce 
the necessary amount of dust mass at f ~ fg„, = 0.6 Gyr The fact 
that the dust content is even higher at earlier times is perhaps 
suggestive of the evolution from a relatively long period of dust 
enshrouded star formation to a phase in which stellar + AGN 
driven winds make the galaxy visible and the QSO shine. Similar 
conclusions can be derived by using the predicted dust-to-gas 
ratio (Fig.fTsTl. 

We now turn our attention to the details of the dust pro- 
duction. In Fig. [14] (upper panel) we also show the cumulative 
mass of dust produced by each single channel (including the dust 
growth) as well as the total mass of dust destroyed. As it hap- 
pens in the solar vicinity (Calura et al., 2008), dust growth and 
destruction dominate at times larger than 0. 1 Gyr The QSO pro- 
duction seems to be relatively unimportant. However, we must 
stress that while SNII production stops and the dust growth is 
strongly suppressed after the galactic wind, the QSO dust pro- 
duction rate should be at its maximum, as shown in Fig. [16] at 
the time of the galactic wind. It is also important to stress that 
the QSO dust should be confined to the very central regions of 
the galaxy, whereas the dust ejected by SNe la and AGBs will 
be more diffuse, therefore the QSO is likely to be the dominant 



source right after the galactic wind in the inner regions of the 
galaxy. 

Although similarly high dust masses have been observed in 
other high redshift QSOs (Beelen et al., 2006), J1148+5251is 
interesting because at its redshift in the favored cosmology, the 
star formation, BH growth and dust production can be repro- 
duced only assuming a very fast evolution. Several other authors 
have attempted to reproduce the observed dust mass: some (e.g. 
Dwek et al., 2007a) claim that the host galaxy is too young to 
have had any AGB contribution to the dust, whereas others (e.g. 
Valiante et al., 2009a) claim that AGB stars are actually a viable 
channel. Moreover, the observed dust mass might be consistent 
with QSO production (see Maiolino et al, 2006). However, it is 
rather difficult to make detailed comparison between different 
models because of their different assumptions. 

For instance, Valiante et al. (2009a) do not consider dust 
growth, have a less efficient destruction (since they do not take 
into account SNe la), but adopt much smaller dust depletion ef- 
ficiencies. Dwek et al. (2007a) assumed that, owing to the young 
(< 0.5 Gyr) inferred age of the host galaxy, only SNII con- 
tribute to the dust production and neglected the role of AGB. 
According to stellar evolution theory, instead, intermediate mass 
stars evolve out of the main sequence after just ~ 30 Myr. 
Therefore, we argue that Dwek (2007)'s assumption is too sim- 
plistic (see also Valiante et al., 2009a) and a more detailed mod- 
els of chemical evolution as the one used in this paper is required 
to assess the AGB contribution. On the other hand, our fiducial 
model needs only 0.6 Gyr to reproduce dust, stellar, gas and QSO 
masses at the redshift (6.4) at which we observe Jl 148+5251. 
Such a time lapse would imply a (maximum) formation redshift 
of 13. This is the value that we adopt for our model, so that 
redshift 6.4 corresponds to tg„. We note that the predicted dust 
mass for our fiducial model crosses the region of the observed 
values in Fig. [14] also at earlier times (0.05 Gyr) which would 
imply - in the current favoured cosmology - a (minimum) for- 
mation redshift of 6.8. Finally, according to Fig. 2, the predicted 
SFR crosses the observed value at ~0.4 Gyr, hence implying a 
formation redshift of ~10. 

In a sense, our predictions give lower limits to the QSO con- 
tribution, which can be higher if either the BH "seed" is larger 
and pre-existing to the galaxy or the growth is super-Eddington 
in the very early stages and quieter near the time of the galactic 
wind. Such a situation is depicted in the lower panel of Fig. [14] 
where the fiducial QSO production (dot-dashed line) adopted 
here is compared to other cases to show the sensitivity to our 
different assumptions. The main properties are summarized in 
Table[T] Case A (solid line) shows the case in which we consider 
the metallicity of the region where the QSO dust forms 0.8 dex 
higher than the mean metallicity predicted by our model. This 
is because BLRs exhibit metallicities up to 10 times the solar 
value, whereas our model predicts the average gas metallicity of 
the central zone to be twice than solar. Such a change in metal- 
licity increases the QSO dust mass by about one order of mag- 
nitude. A similar change can be achieved if we take into account 
the fact that the BH-to-stellar mass ratio is a factor of 2 higher at 
high redshift with respect to the local one (Case B, dotted line, 
e.g. Walter et al., 2004, Maiolino et al., 2007, Lamastra et al., 
2010, Merloni et al., 2010, Wang et al. 2010, but see Willott et 
al., 2010). To take into account the uncertainty in the BH-to- 
stellar mass ratio at high redshift, we make an additional model 
in which this ratio is a factor of 4 higher than the local value 
(Case C, light dot-dashed line). It is worth noting that in both 
Case B and C the final predicted Mbh is still consistent with ob- 
servations, within the uncertainties. Finally, Case D is a model 
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in which we switch off the dust production from stars and as- 
sume that all the dust seen at f ~ tg„ has been produced by the 
QSO. As it can be inferred from comparing the fiducial model to 
cases A, B, C and D, even if the QSO is the only source of dust, 
there will still be degeneracy between the accretion timescale, 
the initial BH seed mass and the metallicity of the clouds where 
the dust forms. The purpose of model M12caseD is, therefore, 
to give just an example and show that a QSO alone might pro- 
duce enough dust if, for instance, we start from a quite high seed 
mass (IO^Mq), while the accretion timescale (0.49 Gyr) is 10 
times longer than the adopted one in the fiducial case. Such a 
situation mirrors the state of the art: the channels for dust pro- 
duction are many, and their relative role cannot be constrained 
by using only the observed dust mass. 

7.1. The QSO abundance ratios 

D'Odorico et al. (2004) determined C, N and ff-element relative 
abundances in the gas surrounding six QSOs (from the NLR) at 
an average redshift of z ~ 2.4, by studying six narrow associated 
absorption systems in Ultraviolet Visual Echelle Spectrograph 
(UVES) high-resolution spectra. They found five systems with a 
metallicity (measured by C/H) consistent with or above the so- 
lar value. They also found a possible correlation of [N/C] with 
[C/H] and an anticorrelation of [Si/C] vs. [C/H]. They suggested 
that to explain such abundance ratios the chemical enrichment in 
the host galaxy had to be very fast, not longer than 1 Gyr This 
conclusion is in agreement with a previous work by Matteucci & 
Padovani (1993) who reproduced the abundaces measured in the 
BLR of QSOs, in particular the oversolar Fe abundance and solar 
O abundance by means of a model of a large elliptical suffering 
a strong burst of star formation. In particular, the high [Fe/Mg] 
ratio claimed in some objects was interpreted as the value in the 
gas around the QSO after star formation has stopped. In fact, in 
such a situation, ff-elements are no more produced whereas Fe 
and N, formed in low and intermediate mass stars, are contin- 
uously produced, thus creating a situation where N and Fe are 
oversolar and the a elements are solar. Our model predicts over- 
solar [Mg/Fe] for the dominant stellar population and underso- 
lar [Mg/Fe] in the gas after the galactic wind, in particular the 
average ratios after the galactic wind in the gas are [Mg/Fe] — 
0.2 dex. That explains also the almost constant metal content 
in QSOs at all redshifts and consequently the similar metallici- 
ties found in BLR at high redshift and in NLR at lower redshift 
(e.g. Nagao et al., 2006a, 2006b, Juarez et al., 2009): most of 
the evolution in these objects has occurred at very ealy epochs 
where SNe II produced most of the metals. After star forma- 
tion stopped due to either galactic winds (Matteucci & Padovani 
1993) or AGN feedback (Romano et al. 2002), the SNe la pro- 
duced the bulk of Fe and low mass stars the bulk of N which 
kept increasing for a while and then settled on a rather constant 
value due to the decrease of Fe an N producers with time (no star 
formation takes place). In Fig. 17 we show our predictions (solid 
line) for [N/C] and [Si/C] versus D'Odorico et al.'s data. There 
is good agreement with the observational trends. A marginal 
agreement (within 3cr) with the high [C/H] observed in QSOs 
(exceeding ~ +1.4 dex) can be obtained only by adopting the 
Maeder (1992) yields with mass loss from massive stars (dot- 
ted lines), which allow the model to reach [C/H] values slightly 
larger than 0.8 dex. With normal yields (e.g. Woosley & Weaver 
1995) we would no go beyond [C/H]=+0.5 dex. 

This result confirms that of Cescutti & al. (2009) who found 
that the Galactic bulge high C abundances can be reproduced 
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Fig. 14. Upper panel: Dust mass (Red solid line) as function 
of galactic age. The cumulative contribution to the dust mass by 
SNIa, SNII, AGB and QSO (fiducial case)) as well as the integral 
of the the dust growth and destruction are shown. The observed 
values are given by the shaded area. The arrow shows the time 
at which we predict the conditions for the galactic wind to occur 
in most of the galaxy. The model refers to a galaxy of IO'^Mq 
(model M12d). Lower panel: as above, in this case only the total 
dust and the QSO contribution are plotted. The QSO fiducial 
dust production is compared to a case in which we assume the 
metallicity of the dust forming regions 0.8 dex higher than the 
average (Case A), and to two cases in which the relation between 
BH mass and stellar mass is "super-Magorrian" (Cases B and C, 
see text). Finally, we present a model (Case D) in which the QSO 
is fine-tuned to be the only source of dust. 

only by adopting yields from massive stars with mass loss 
strongly dependent on metallicity. 

8. Conclusions 

In this paper we have studied the evolution of several chemical 
elements (C, N, O, Si, Mg, Fe) both in the gas and in the dust in 
elliptical galaxies of different masses (from lO'^Mo to IO'-Mq). 
Our aim was to study the evolution of gas and dust at high red- 
shift and the effect of different environments. We assumed that 
dust is produced in: i) SNe II, ii) SNe la, iii) AGB stars and iv) 
QSOs. We took into account also dust destruction and accretion. 
These models will help up-coming observations of distant ob- 
jects in quantifying their morphologies, past SFR and mass by 
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Fig. 15. The predicted dust to gas mass ratio as a function of 
galactic age for model Ml 2d (continuous line) compared to the 
observational range of data (shaded area). 




0.2 0.4 0.6 0.2 0.4 0.6 

galactic age (Gyr) galactic age (Gyr) 



Fig. 16. Dust production rate of SNe la (left upper panel), SNe 
II(left lower panel), AGB stars (right upper panel), QSO(right 
lower panel) as functions of galactic age for different elements 
as indicated in the figure, for model Ml 2d. 

means of abundances and abundance ratios, as well as the role 
of the dust in the measurement of the above quantities. In partic- 
ular, our general predictions are that: 

- The downsizing implies that at [Fe/H]< -1 the most mas- 
sive objects have lower [N/Fe] ratios than the lower mass 
ones, irrespective of the presence of dust. The same is true 
for the [N/O] ratios in a broader range of metallicities. 

- A dust mass-stellar mass relation exists, with more mas- 
sive galaxies attaining a higher dust content at earlier time. 
Thererefore the expected evolution is similar to the observed 
high redshift evolution of the mass-metallicity relation. 

These predictions await confirmation from future observations. 
In order to test the model we compared our results for the lower 
mass galaxies with LBG galaxies and for the highest mass limit 
with QSOs. In particular, we compared our model for LBGs also 
with the specific galaxy MS 1512-cB58. Unfortunately, the ob- 
served abundance ratios from diff'erent authors show very dif- 
ferent values. Our predicted [N/Fe] ratio assuming no primary 
production of N from massive stars and with dust reproduces 



Fig. 17. Predicted and observed [N/C] and [Si/C] ratios versus 
[C/H]. The model refer to an elliptical of IO'^Mq. In particular, 
the solid line refers to the fiducial case (model Ml 2d), whereas 
the dotted line gives the predictions in which the metallicity de- 
pendent yield from Woosley & Weaver (1995) + Maeder (1992) 
are adopted (model M12MDd). Data: NLR by D'Odorico et al. 
(2004). 



very well tha data of Pettini & al. (2002). On the other hand, the 
model assuming dust and primary N production from massive 
stars predicts the [N/Fe] ratio in very good agreement with the 
data of Teplitz et al. (2000). Therefore, no conclusions can be 
drawn on this point. 

As far as these specific class of objects are concerned, our 
findings can be summarized as follows: 

- LBG galaxies are likely to be young ellipticals of interme- 
diate mass(10"' - 3 ■ lO^^M©) experiencing moderate SFR 
and galactic winds. Our elliptical model for 3 ■ IO^'Mq well 
reproduces the [O/H] abundance as a function of redshift for 
the LBG galaxies. Our model implies a redshift of formation 
for these galaxies in the z - 2.3 - 4 range. 

- The predicted [Mg/Fe] and [Si/Fe] ratios have been com- 
pared with data for MS 1512-cB58 by Pettini & al. (2002) 
and the theoretical values are lower than the observed ones, 
perhaps indicating that too low a fraction of Fe is predicted 
to be in dust. 

- From the comparison of our model results for a massive el- 
liptical (M = lO'-Mo) with one of the most distant QSO ever 
observed J 1 148+5251, we derived a redshift of formation for 
this object z > 7 (and likely larger than 10). 

- The predicted total amount of dust is in good agreement with 
the dust seen in this QSO. We tested also the hypothesis 
that the QSO itself produced dust but this production appears 
negligible compared to that from stellar sources, unless one 
focuses on the very central regions at times very close to the 
galactic wind onset. Future observations can provide better 
constraints to the QSO dust model. 

- We also compared our model results for QSO hosts with data 
from NLR in QSO hosts and we found a good agreement for 
[N/C] versus [C/H] and [Si/C] versus [C/H]. The very high 
C abundance observed in these QSOs can be explained only 
by assuming yields with mass loss from massive stars with a 
strong dependence on metallicity, as those of Maeder (1992). 
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